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We study screening of optical singularities in random optical fields with two widely different length 
scales. We call the speckle patterns generated by such fields speckled speckle, because the major 
speckle spots in the pattern are themselves highly speckled. We study combinations of fields whose 
components exhibit short- and long-range correlations, and find urmsual forms of screening. 
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I. INTRODUCTION 

Screening of charged topological singularities - vortices 
[1, 2 (Chap. 5), 3 (Sect. 4.8)] in scalar fields, C points 
[2, (Chaps. 12 and 13)] in vector fields - has been exten- 
sively studied in random fields with a single correlation 
length [4—19]; here we study screening of these singulari- 
ties in random fields with two widely different correlation 
lengths. We call such fields "speckled speckle" because, 
as illustrated in Fig. IT|a), the major speckle spots of the 
field are themselves highly speckled. Speckled speckle 
fields can be generated by illuminating a random diffuser 
with two concentric, overlapping beams: one, the a beam, 
is tightly focused and intense, the other, the b beam, is 
weak and diffuse. 

The statistical properties of speckled speckle can be 
highly anomalous, with relative number densities of crit- 
ical points (vortices, C points, extrema, and umbilic 
points) differing from normal speckle values by orders 
of magnitude [20,21]. The spatial arrangement of vor- 
tices and C points is also anomalous, with these singular- 
ities forming dense clusters of a kind not found in normal 
speckle fields. Fig. [ija) [20,21]. 

Screening can be either short- or long-ranged. Nonsin- 
gular random sources produce random fields that exhibit 
short range screening [4—19]. In such systems positive 
(negative) topological charges are surrounded by a local 
net excess of negative (positive) charge, leading to charge 
neutrality within a characteristic distance, the screening 
length, that can be less than the average separation be- 
tween charges [19]. 

Singular sources, such as a ring of finite radius but 
zero width [8] , produce random fields that exhibits long- 
range screening [8, 16, 19]. The singularities in the field 
produced by a ring form a quasi-lattice in which pos- 
itive/negative singularities occupy alternate corners of 
a square cell. Fig. fljb). Local defects in the lattice 
destroy the local charge neutrality that would produce 
short range screening, and screening sets in only asymp- 
totically. 

Thus, both short- and long-range screening depend 
upon the spatial arrangements of the charges. These 
arrangements are anomalous in speckled speckle, which 
can therefore be expected to exhibit unusual forms of 
screening. 

A major observable consequence of screening is strong 
damping of fluctuations of the topological charge Q. 



These fluctuations are characterized by their variance 
(Q^); the behavior of (Q^) for speckled speckle is our 
main concern here. 

The plan of this paper is as follows. In Section II we 
discuss the charge variance in a bounded region, in Sec- 
tion III we review {Q^} in normal speckle for fields with 
short- and with long-range correlations, and in Sections 
IV- VII we present results for (Q^) for four qualitatively 
different forms of speckled speckle. We briefly summa- 
rize our main findings in the concluding Section VIII. 
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FIG. 1: Vortex structures. Positive (negative) vortices are 
shown by white (black) filled circles, (a) Speckled speckle. A 
random diffuser is illuminated by two concentric disks of light, 
a and b. The diameter of disk b is ten times the diameter of 
a; the total optical power in a, however, is 10 times that in 
b. Major (minor) speckle spots in the speckled speckle field 
are due primarily to beam a (6). Vortices of the combined 
beam cluster in the dark regions between a field speckle spots 
because they require perfect destructive interference between 
the strong a and weak 6 fields. In normal speckle produced 
by a single disk vortices tend to be uniformly distributed with 
only minor clustering, (b) Normal speckle phase map pro- 
duced by a random diffuser illuminated with a single ring; the 
vortices tend to form a lattice with a square unit cell. 



II. CHARGE VARIANCE IN SHORT- AND 
LONG-RANGE SCREENING 



Source Distributions and Autocorrelation 
Functions 



We assume isotropy, circular Gaussian statistics [3 
(Chap 2), 22 (Chap. 2)], and a circular region of radius 
R. The charge variance (Q^) in this region is related 
to the autocorrelation function of the field W{r) by [19 
(Eq. 39)], 
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where W'{r) — dW{r)/dr. The number density of 
charges i] is [23] 
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For the short range correlations produced by extended 
sources, W'{r) decays rapidly with r, and in the limit of 
large R, 
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a 2.81182, 


(4b) 
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a 0.563047, 


(4c) 



Thus, for short-range screening (Q'^) grows with the 
perimeter, i.e. (Q^) ^ R ^ \N, in contrast to the 
case of no screening, where (Q^) grows with the area, 
{Q^)^R^r^ N [8,19]. 

For the long-range screening produced by a singular 
ring source of radius p and zero width, the large R limit 
is [19 (Eq. 48)], 






V = 



where 7 « 0.577216 is Euler's constant. Thus, long- 
range screening yields a charge variance that grows 
asymptotically as R In R; this growth rate is significantly 
faster than the short-range growth rate proportional to 
R, but is very much slower than the unscreened rate pro- 
portional to R^. For, say, a large region that contains 
10^ charges, short-range (long-range) screening damps 
out charge fluctuations relative to no screening by a fac- 
tor of - 100 (-22). 



III. CHARGE VARIANCE IN NORMAL 
SPECKLE 

We review here the charge variance in normal speckle 
produced by sources with a single characteristic length 
scale. In later sections we build our composite sources 
with their two different length scales from binary combi- 
nations of these single sources, and compare composite- 
source charge variances with single-source variances. 



Listed below are the source distributions S{u), where 
u measures radial displacements in the source plane, the 
total optical power in each source, P, and the autocor- 
relation functions W{r) of the speckle field. S{u) and 
W{r) are related by the VanCittert-Zernike theorem [3 
(Sect. 4), 22 (Sect. 5.6)]. We study four fields with 
autocorrelation functions that decay at different rates. 

(i) A Gaussian, superscript (G), of 1/e width p, and in- 
tensity (optical power/unit area) I^'^' at the peak; hence- 
forth the "Gaussian" , 



^(G)(u)^/(G)exp(_[M/(2p)]' 

p(G)=4V/^°', 
W^'^^r) =exp(-p2r2) . 



(5a) 

(5b) 
(5c) 



(ii) A uniform disk, superscript (D), of radius p, and 
uniform intensity J'^'^-' ; henceforth the "Disk" , 



5(D)(u)=/(D)e(w-p), 

W'-^\r) = 2Ji{pr)/{pr). 
function 



(6a) 
(6b) 
(6c) 
defined by 



O {x) is the Heaviside step 
e(j:<0) = 0, e(a;>0) = 1. 

{in) A nonuniform (inverse square root) disk, super- 
script (S), of radius p, with intensity J'^^-' at the disk 
center, 

/-(s) 
^(S) (u) = 9 {u~p), (7a) 

VI - {u/pf 
F(S)^27r//(S), (7b) 

W^(s)(r) = /(S)sinc(pr), (7c) 



where, sinc(a;) = sin (x) jx. In what follows we refer to 
this source as the "Sine" . 

(iv) A singular ring, superscript (R), of radius p, which 
we write as 



p(R) 

l^(R)(r) 



/(«)e,5(w- 

27rpe/(^), 

Mpr), 



P) 



(8a) 
(8b) 
(8c) 



where 5 {x) is the Dirac delta function. In what follows 
we refer to this source as the "Ring" . 

S'^^^u) is the limit of a finite width annulus s^^^ (u) 
of mean radius p and width e. 



s{u) = e{u~p~e/2) 
lim [s (u) /e] 



e{u-p + e/2), (9a) 

:6{u~p). (9b) 



/' ^ is therefore the uniform intensity in th e an nulus. In 
the limit e ^ 0, /^^' diverges, p(^) in Eq. (8b), however, 
is assumed to remain finite. 



W and (Q^) for the finite width annulus {e > 0) is 
discussed in [19 (Sect. 5)], where it is shown that over 
the region er < I the experimentally attainable annulus 
is an excellent approximation to the theoretical singular 
ring. 



B. Charge Variance 

In Fig. 0we plot (Q^) / (pR) vs. pR for the above four 
sources. Tor the Gaussian (G), Disk (D), and Sine (S), 
screening is short-ranged, and the large R limit of (Q^) 
is given in Eq. dsl. 

For the Gaussian, Eq. ^ can be evaluated analyti- 
cally, yielding 

(Q^) Gauss ~ \ ^2AC (3/2) pR = 0.521093pi?, (10) 

with Q [x) is the Riemann zeta function, whereas for the 
Disk (D) and Sine (S), Eq. (IS]) is evaluated numerically, 
yielding 



0.6 



{Q') 



Disk 



and 



{Q\: 



Si nc 



0.227210pi?, 



0.305898pi?. 



(11) 



(12) 



For the Ring (R) screening is long-ranged, and the 
large R limit of (Q^) / {pR) is given in Eq. H. 
For very small R for all sources [19 (Eq. 82j], 



{Q') 



R^O 



TjTrR^ 



N. 



(13) 



i.e. there is no screening; the reason is that for a suffi- 
ciently small area the probability of finding the required 
screening charges within the area is vanishingly small. 
This result is illustrated in Fig. |2jc) for all four sources. 



IV. COMPOSITE SOURCES AND THEIR 
AUTOCORRELATION FUNCTIONS 

For scalar (single component) fields the a and b beams 
have the same, say, linear polarization, and the singu- 
larities whose screening is of interest here are the phase 
vortices [1 — 19]. For vector (two component) fields the 
a and b beams have orthogonal linear polarizations, and 
the relevant singularities that screen each other are either 
right- or left-handed C points [2] : right-handed C points 
do not screen left-handed ones, and vice versa. 

We write our composite source as 



5(T"^^) (.) = ^f) (.) + sr (-) , 



(14) 



where the source type specifier (T^Tft) is a binary combi- 
nation of Ta,6 = G, D ,S, R. From the VanCittert-Zernike 
theorem, the corresponding autocorrelation function is 



M/(T.T,)(^) 



WP (r) + 



i^(T„T,)^m (^) 



1 + K(TaT.) 



(15) 
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FIG. 2: Charge variance (Q^) for normal speckle obtained 
from numerical integration of Eq. ([l|. (a) Short-range 
screening, (b) Long-range screening. The dependence on 
the parameter p in Eqs. ([5|-([8| is here scaled out by plot- 
ting (Cf'^J / {pR) vs. pR. As can be seen, for pR > 1, the 
result s in ( a) quickly asymptote to the theoretical values in 
Eqs. ( |10[ |-( 12 I, whereas the result in (b) asymptotes to the 
theoretical form (thin straight line) in Eq. |4|. (c) (Q^) for 
small R. In all four cases (G - Gauss, D - Disk, S - Sine, and 
R - Ring), the curves approach the i? — > limit given in Eq. 
(131. 



where the dimcnsionlcss constant 

^(T„T,)^pm/pm. (16) 

sSi"")' ^i?' ^^d ^i? ('^)' ^^e listed ™ Eqs. Q-^, 
with p = a,b as appropriate. 

Similarly, the number density of singularities is for 
composite scalar fields, 

2tT 1 + if (TaT,) 



rj'- 



(17) 



This is also the number density of right- and of left- 
handed C points. 



As a specific example, for the intense, tightly focused 
a beam a Gaussian (G), and the diffuse weak b beam a 
Sine (S), 



Sl''\u) = lPe^p[-[u/i2a)f), 



AS) 






^(GS) ^ 



exp (— ar) + iiT^'^^^sinc (br) 



1 + if (GS) 



2 AS) 



2VL 



I 



a'lP^ 



V 



(CS) _ 2a^ + Ki^^^by3 



1 + if (GS) 



(18a) 
(18b) 

(18c) 
(18d) 

(18e) 



In the composite-source examples that follow we usu- 
ally take K ^ 0.01 and b/a = 100. There are two reasons 
for these choices: (i) small K together with large b/a pro- 
duces results that vividly illustrate the unusual screening 
properties of speckled speckle, and (ii) this combination 
of parameters permits a significant degree of analysis. 
We consider that these parameters, which are convenient 
for the theoretician, to be experimentally possible; ad- 
mittedly, they may be difficult to achieve in practice. 

With the above choice of parameters, the R depen- 
dence of \Q'^) separates into three distinct regions: 

I. R< 1/5. 



In this region Eq. ( 13 1 holds for all composite sources 



with rj equal to the number density of b field charges 



K 



V-Vb 



'2^^^^ 



"(T) 



(0). 



(19) 



The reason is that for sufficiently small R the probabil- 
ity of finding an a beam charge in the area is negligible; 
only b beam charges are present, so only these charges 
contribute to (Q^). This result is verified by direct com- 
parison (not shown) with the exact result in Eq. (nj). 

II. l/b<R< 1/a. 

In this region, a in Eq. (151 may be set equal to zero, 
because for ar < 1, Wa{r) w Wa{0) = 1, and W^{r) « 
W^{0) = 0, independent of the a beam type - G, D, S, or 
R. As will become apparent, these approximations yield 
good agreement with the exact result in Eq. (fl|. 

III. R > 1/a. 

In this region there is no accurate approximation that 
is applicable, however, as discussed below, approxima- 
tions good to ~ 10% are available. 

Below we discuss screening in composite beams for the 
four qualitatively different combinations of a and b beams 
in which the individual beams exhibit either short- or 
long-range correlations. 



V. BOTH BEAM a AND BEAM b EXHIBIT 
SHORT-RANGE SCREENING 

We start with region II, 1/6 < R < 1/a. Here the area 
contains many b charges but practically no a charges. We 
denote the b charge contribution to (Q^) by {Q'^)i^, where 



{Q'\ 



R_ 

2^' 



if(T„T. 



dw'i^\r)ldr 



l-M/f)(r) 



-dr. (20) 



In obtaining this result we make use of the fact that K <C 

2 



1 and 



KW, 



(T) 



<M^, 



(T) 



For b a Gaussian (G), Eq. (20 1 can be evaluated ana- 
lytically, and we have. 



m 



(G) 
h 



2^ 
0.454843ii:('^»°)5i?. 



i?, 



(21) 



cally, and we have for the Disk (D) 



For b a Disk or a Sine, Eq. (|20|) is evaluated numeri- 

(22) 



1\{D) 



and for the Sine (S), 



0.187153is:('^"°)6i?, 



0.238531if(^"^)&i?. 



(23) 



In region III, R > 1/a, we assume that the b charges 
continue to contribute (Q^)^ to (Q^). In addition, the 
area now includes many a charges. We label the contri- 
bution of these charges (Q^) , and write 



{Q')-{Q')a + {Q\ 



(24) 



i.e., we neglect interactions between the a and b beams, 
and the possibility that a charges can screen b charges, 
and vice versa. Within the framework of this approxi- 
mation we write. 



.2\(G) 



(Qyr = {Q')aa. 



[Eq.dOl] 



-,2\(-D) 



{Qy:' = {Q')n.s.' [Eq-©] 



.2\(S) 



(QX' = (Q') 



Sine ' 



Eq.(12l] 



(25) 



(26) 



(27) 



with p replaced by a. 

We illustrate the above in Fig. |3]. As can be seen, Eqs. 



(21l-(23) are good approximations to the exact results. 



whereas Eqs. (24)-(27l are good only to order of 10%. 



We note that the well defined steps in this figure pro- 
vide striking visual confirmation of the fact that there are 
two widely different length scales. The first step starts, 
as expected, at i? ~ 1/b, the second at i? ~ 1/a. 
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FIG. 3: (Q^) jR vs. 7? for composite sources in wiiicii both 
the a and h beams have short range correlations. The thick 
curves are the exact result in EqJTJ . The solid lines labeled 
B are the theory in Eqs. (20 1- ( |23[), the dashed lines labeled 
A are the theory in Eqs. (24 1- ( |27| , whereas the solid lines 
labeled C are the short range screening result, Eq. ([3|. (a) 
GG, beams a and 6 are Gaussians. Beam parameters are 
a = 1,6 = 100,74' = 0.01. (b) SG, beam a is a Sine, beam 
6 is a Gaussian. Beam parameters are the same as in (a), 
(c) DG, beam a is a Disk, beam & is a Gaussian. Beam 
parameters are a= 1, & = 10, K = 0.04. In all three examples 
displayed here the h field was chosen to be a Gaussian in order 
to emphasize that the agreement with the theory represented 
by line B does not depend on the nature of the o beam. Other 
short range choices for the h beam. Disk or Sine, show equally 
good agreement. 



Neglecting again cross screening of a and h charges, the 
total charge variance in region III is approximated by the 
sum of a and h beam contributions, Eq. (24). 

In Fig. [ijwe plot (Q^) jji vs. R for RG and DG, again 
obtaining m the region of short range screening, region 
B, a plateau that is in good agreement with the calcu- 
lated value for {Q'^\. As expected, in region A, the 
region of long range screening, (Q^) /R grows linearly 
with In R. We note that here not only is the calculated 
slope, Eq. (|4|, in close agreement with the exact result 
(thick curves), but also that Eq. ( p4| ) provides a rather 
reasonable description of the data. Similar good agree- 
ment is obtained for RS (not shown). 



VII. BEAM a EXHIBITS SHORT-RANGE 

SCREENING, BEAM b EXHIBITS LONG-RANGE 

SCREENING 



In region II, 1/6 < R < 1/a, we again set a = in 
pj/(TaT6)^ and obtain for beam b a Ring, 
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VI. BEAM a EXHIBITS LONG-RANGE 

SCREENING, BEAM h EXHIBITS 

SHORT-RANGE SCREENING 



As discussed in the previous section, in region II, 1/6 < 



R < 1/a, (Q^) is dominated by (Q^)^, Eqs. (21 1- (23 1 



In region III, where R > 1/a, beam a exhibits long-range 
screening, and (Q^) is given by Eq. mh with p = a. 



FIG. 4: (Q^) /R vs. R for composite sources in which the 
a beam has long-range and the b beam short-range correla- 
tions. The thick curves are the exact result in Eq. ([T]). The 
solid line s labe led B are the short-range screening theory in 
Eqs. ( 20 1- ( 22 I , the dashed lines labeled A are the long-range 
screening result in Eq. (H. (a) RG, beam a is a Ring, 6 a 
Gaussian, (b) RD, beam a is a Ring, beam b a Disk. Beam 
parameters in both (a) and (b) are a — l,b = 100, K — 0.01. 



As before, we have assumed K ^ 1, 



KW, 



(T) 



<$:w, 



(T) 



Proceeding as in [19, (Eqs. 40-48)], we obtain 



(Q') 



KhR 



[T + ln(6i?)] , 



b 27r2 
JP = 7rX+7 + 51n2-3« 4.84258, 

Jo {x) Jl {x ) 
1 - Jo {x) 



-dx « 1.20946. 



(29a) 
(29b) 

(29c) 



We illustrate the above in Fig. [sj In region B, 10 ^ < 
R < 10-\ for small K (0.015) Eqs. J29]) provide a good 
description of the exact result. Fig. Isfa), whereas when 
K is no longer small (0.2) the expected deviations appear. 
Fig. gb). 

In Fig. [sja) a plateau appears for i? > 1, apparently 
consistent with the fact that short-range screening a field 
charges become important in this region. But this appar- 
ent plateau is misleading, because the long-range screen- 
ing of the h field charges can never saturate: regardless of 
the nature of the a field, (Q^) /R for the h charges must 
diverge logarithmically for large R. The rate (slope) 
of this logarithmic divergence, however, is K dependent, 
being small. Fig. Isja) (large. Fig. Isl^b)) for small (large) 
K. 



VIII. BOTH BEAM a AND BEAM b EXHIBIT 
LONG-RANGE SCREENING 

The case of two Rings is illustrated in Fig. [6] As 
before, Eq. (29 1 holds for the b field charges. But now 



the a field charges also exhibit long-range screening, and 
(Q^) /R exhibits the expected large R logarithmic diver- 
gence also for small K (0.015). 




FIG. 6: (Q^) / R vs. R for a composite source in which both 
beam a and beam b are Rings. Beam parameters are a = 
1,6= 100, K — 0.015. The thick curve is the exact result in 
Eq. (11]). The solid hne labeled B is the long-range screening 
theory in Eqs. pgl. 




FIG. 5: (Q^) /R vs. R for composite sources in which beam 
a is a Gaussian with short-range correlations, and beam b 
a Ring with long-range correlations. Beam parameters are 
a = l,b = 100. The thick curves are the exact result in Eq. 
([l]). The solid lines labeled B are the long-range screening 
theory in Eqs. (|29|. (a) K = 0.015. (b) K = 0.2. 



IX. SUMMARY 

Screening in speckled speckle produced by composite 
sources has been studied in the limit of two widely dif- 
ferent characteristic length scales in the speckle pattern 
and two widely different intensities in the sources. The 
interesting combination, emphasized here, is an intense 
field with a long length scale (the a field) perturbed by 
a weak field that has a short length scale (the b field) . 

Exact, and approximate, results have been presented 
for all four combinations of short- and long-range screen- 
ing. When the correlations in both fields are short-range, 
cross screening between a and b field singularities appears 
to be of only secondary importance, and the results of the 
exact calculation are found to decompose, approximately, 
into a sum of screening contributions, one for each field. 
Fig. [3j A similar decomposition is found to hold if the 
b field exhibits short-range screening and the a field ex- 
hibits long-range screening. Fig. |4J If the b field exhibits 
long-range screening, however, then no decomposition is 
vahd. Figs. [5] and |6] 

The unusual screening properties of speckled speckle 
are most pronounced when the ratio (b field to a field) of 
length scales is 100 : 1 or greater, and the ratio of optical 
powers in the sources (the parameter K) is 1 : 100 or 
less. Less extreme parameter ratios can also yield useful 
results, however, as in Fig. ^b), so that the experimen- 
tal study of screening in speckled speckle appears to be 
entirely feasible. 
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